Pathogenic TNNI1 variants disrupt sarcomere contractility resulting in hypo- and hypercontractile muscle disease.

Troponin I (TnI) regulates thin filament activation and muscle contraction. Two isoforms, TnI-fast (TNNI2) and TnI-slow (TNNI1), are predominantly expressed in fast- and slow-twitch myofibers, respectively. TNNI2 variants are a rare cause of arthrogryposis, whereas TNNI1 variants have not been conclusively established to cause skeletal myopathy. We identified recessive loss-of-function TNNI1 variants as well as dominant gain-of-function TNNI1 variants as a cause of muscle disease, each with distinct physiological consequences and disease mechanisms. We identified three families with biallelic TNNI1 variants (F1: p.R14H/c.190-9G>A, F2 and F3: homozygous p.R14C), resulting in loss of function, manifesting with early-onset progressive muscle weakness and rod formation on histology. We also identified two families with a dominantly acting heterozygous TNNI1 variant (F4: p.R174Q and F5: p.K176del), resulting in gain of function, manifesting with muscle cramping, myalgias, and rod formation in F5. In zebrafish, TnI proteins with either of the missense variants (p.R14H; p.R174Q) incorporated into thin filaments. Molecular dynamics simulations suggested that the loss-of-function p.R14H variant decouples TnI from TnC, which was supported by functional studies showing a reduced force response of sarcomeres to submaximal [Ca2+] in patient myofibers. This contractile deficit could be reversed by a slow skeletal muscle troponin activator. In contrast, patient myofibers with the gain-of-function p.R174Q variant showed an increased force to submaximal [Ca2+], which was reversed by the small-molecule drug mavacamten. Our findings demonstrated that TNNI1 variants can cause muscle disease with variant-specific pathomechanisms, manifesting as either a hypo- or a hypercontractile phenotype, suggesting rational therapeutic strategies for each mechanism.

Peptide-scFv antigen recognition domains effectively confer CAR T cell multiantigen specificity.

The emergence of immune escape is a significant roadblock to developing effective chimeric antigen receptor (CAR) T cell therapies against hematological malignancies, including acute myeloid leukemia (AML). Here, we demonstrate feasibility of targeting two antigens simultaneously by combining a GRP78-specific peptide antigen recognition domain with a CD123-specific scFv to generate a peptide-scFv bispecific antigen recognition domain (78.123). To achieve this, we test linkers with varying length and flexibility and perform immunophenotypic and functional characterization. We demonstrate that bispecific CAR T cells successfully recognize and kill tumor cells that express GRP78, CD123, or both antigens and have improved antitumor activity compared to their monospecific counterparts when both antigens are expressed. Protein structure prediction suggests that linker length and compactness influence the functionality of the generated bispecific CARs. Thus, we present a bispecific CAR design strategy to prevent immune escape in AML that can be extended to other peptide-scFv combinations.

Recent breakthroughs in computational structural biology harnessing the power of sequences and structures.

Recent advances in computational approaches and their integration into structural biology enable tackling increasingly complex questions. Here, we discuss several key areas, highlighting breakthroughs and remaining challenges. Theoretical modeling has provided tools to accurately predict and design protein structures on a scale currently difficult to achieve using experimental approaches. Molecular Dynamics simulations have become faster and more precise, delivering actionable information inaccessible by current experimental methods. Virtual screening workflows allow a high-throughput approach to discover ligands that bind and modulate protein function, while Machine Learning methods enable the design of proteins with new functionalities. Integrative structural biology combines several of these approaches, pushing the frontiers of structural and functional characterization to ever larger systems, advancing towards a complete understanding of the living cell. These breakthroughs will accelerate and significantly impact diverse areas of science.

Prediction of Kv11.1 potassium channel PAS-domain variants trafficking via machine learning.

Congenital long QT syndrome (LQTS) is characterized by a prolonged QT-interval on an electrocardiogram (ECG). An abnormal prolongation in the QT-interval increases the risk for fatal arrhythmias. Genetic variants in several different cardiac ion channel genes, including KCNH2, are known to cause LQTS. Here, we evaluated whether structure-based molecular dynamics (MD) simulations and machine learning (ML) could improve the identification of missense variants in LQTS-linked genes. To do this, we investigated KCNH2 missense variants in the Kv11.1 channel protein shown to have wild type (WT) like or class II (trafficking-deficient) phenotypes in vitro. We focused on KCNH2 missense variants that disrupt normal Kv11.1 channel protein trafficking, as it is the most common phenotype for LQTS-associated variants. Specifically, we used computational techniques to correlate structural and dynamic changes in the Kv11.1 channel protein PAS domain (PASD) with Kv11.1 channel protein trafficking phenotypes. These simulations unveiled several molecular features, including the numbers of hydrating waters and hydrogen bonding pairs, as well as folding free energy scores, that are predictive of trafficking. We then used statistical and machine learning (ML) (Decision tree (DT), Random forest (RF), and Support vector machine (SVM)) techniques to classify variants using these simulation-derived features. Together with bioinformatics data, such as sequence conservation and folding energies, we were able to predict with reasonable accuracy (≈75%) which KCNH2 variants do not traffic normally. We conclude that structure-based simulations of KCNH2 variants localized to the Kv11.1 channel PASD led to an improvement in classification accuracy. Therefore, this approach should be considered to complement the classification of variant of unknown significance (VUS) in the Kv11.1 channel PASD.

An investigation of the YidC-mediated membrane insertion of Pf3 coat protein using molecular dynamics simulations.

YidC is a membrane protein that facilitates the insertion of newly synthesized proteins into lipid membranes. Through YidC, proteins are inserted into the lipid bilayer via the SecYEG-dependent complex. Additionally, YidC functions as a chaperone in protein folding processes. Several studies have provided evidence of its independent insertion mechanism. However, the mechanistic details of the YidC SecY-independent protein insertion mechanism remain elusive at the molecular level. This study elucidates the insertion mechanism of YidC at an atomic level through a combination of equilibrium and non-equilibrium molecular dynamics (MD) simulations. Different docking models of YidC-Pf3 in the lipid bilayer were built in this study to better understand the insertion mechanism. To conduct a complete investigation of the conformational difference between the two docking models developed, we used classical molecular dynamics simulations supplemented with a non-equilibrium technique. Our findings indicate that the YidC transmembrane (TM) groove is essential for this high-affinity interaction and that the hydrophilic nature of the YidC groove plays an important role in protein transport across the cytoplasmic membrane bilayer to the periplasmic side. At different stages of the insertion process, conformational changes in YidC's TM domain and membrane core have a mechanistic effect on the Pf3 coat protein. Furthermore, during the insertion phase, the hydration and dehydration of the YidC's hydrophilic groove are critical. These results demonstrate that Pf3 coat protein interactions with the membrane and YidC vary in different conformational states during the insertion process. Finally, this extensive study directly confirms that YidC functions as an independent insertase.

Elucidating the molecular basis of spontaneous activation in an engineered mechanosensitive channel.

Mechanosensitive channel of large conductance (MscL) detects and responds to changes in the pressure profile of cellular membranes and transduces the mechanical energy into electrical and/or chemical signals. MscL can be activated using ultrasonic or chemical activation methods to improve the absorption of medicines and bioactive compounds into cells. However, re-engineering chemical signals such as pH change can trigger channel activation in MscL. This study elucidates the activation mechanism of an engineered MscL at an atomic level through a combination of equilibrium and non-equilibrium (NE) molecular dynamics (MD) simulations. Comparing the wild-type (WT) and engineered MscL activation processes suggests that the two systems are likely associated with different active states and different transition pathways. These findings indicate that (1) periplasmic loops play a key role in the activation process of MscL, (2) the loss of various backbone-backbone hydrogen bonds and salt bridge interactions in the engineered MscL channel causes the spontaneous opening of the channel, and (3) the most significant interactions lost during the activation process are between the transmembrane helices 1 and 2 in engineered MscL channel. The orientation-based biasing approach for producing and optimizing an open MscL model used in this work is a promising way to characterize unknown protein functional states and investigate the activation processes in ion channels and transmembrane proteins in general. This work paves the way for a computational framework for engineering more efficient pH-sensing mechanosensitive channels.

Modulation of P2X4 pore closure by magnesium, potassium, and ATP.

The P2X4 receptor plays a prominent role in cellular responses to extracellular ATP. Through classical all-atom molecular dynamics (MD) simulations totaling 24 µs we have investigated how metal-complexed ATP stabilizes the channel's open state and prevents its closing. We have identified two metal-binding sites, Mg2+ and potassium K+, one at the intersection of the three subunits in the ectodomain (MBS1) and the second one near the ATP-binding site (MBS2), similar to those characterized in Gulf Coast P2X. Our data indicate that when Mg2+ and K+ ions are complexed with ATP, the channel is locked into an open state. Interestingly, irrespective of the number of bound ATP molecules, Mg2+ ions bound to the MBS2 impeded the collapse of the open-state protein to a closed state by stabilizing the ATP-protein interactions. However, when Mg2+ in the MBS2 was replaced with K+ ions, as might be expected when in equilibrium with an extracellular solution, the interactions between the subunits were weakened and the pore collapsed. This collapse was apparent when fewer than two ATPs were bound to MBS2 in the presence of K+. Therefore, the different capacities of common cations to stabilize the channel may underlie a mechanism governing P2X4 channel gating in physiological systems. This study therefore provides structural insights into the differential modulation of ATP activation of P2X4 by Mg2+ and K+.

Mechanistic Picture for Chemomechanical Coupling in a Bacterial Proton-Coupled Oligopeptide Transporter from Streptococcus Thermophilus.

Proton-coupled oligopeptide transporters (POTs) use the proton electrochemical gradient to transport peptides across the cell membrane. Despite the significant biological and biomedical relevance of these proteins, a detailed mechanistic picture for chemomechanical couplings involved in substrate/proton transport and protein structural changes is missing. Therefore, we performed microsecond-level molecular dynamics simulations of bacterial POT PepTSt, which shares ∼80% sequence identity with the human POT, PepT1, in the substrate-binding region. Three different conformational states of PepTSt were simulated, including (i) occluded, apo, (ii) inward-facing, apo, and (iii) inward-facingoccluded, Leu-Ala bound. We propose that the interaction of R33 with E299 and E300 acts as a conformational switch (i.e., to trigger the conformational change from an inward- to outward-facing state) in the substrate transport. Additionally, we propose that E299 and E400 disengage from interacting with the substrate either through protonation or through coordination with a cation for the substrate to get transported. This study provides clues to understand the chemomechanical couplings in POTs and paves the way to decipher the molecular-level underpinnings of the structure-function relationship in this important family of transporters.

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins.

Members of the parvalbumin (PV) family of calcium (Ca2+) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca2+ affinities and selectivities against competing ions like Mg2+ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca2+ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca2+ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg2+ to these isoforms, which is much weaker than Ca2+. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg2+ from Ca2+, they did not indicate that Ca2+ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca2+ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg2+ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, ß-sheet character, and hydrogen bonding. Hence, Ca2+ affinity and selectivity against Mg2+ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.

Pathogenic variants in TNNC2 cause congenital myopathy due to an impaired force response to calcium.

Troponin C (TnC) is a critical regulator of skeletal muscle contraction; it binds Ca2+ to activate muscle contraction. Surprisingly, the gene encoding fast skeletal TnC (TNNC2) has not yet been implicated in muscle disease. Here, we report 2 families with pathogenic variants in TNNC2. Patients present with a distinct, dominantly inherited congenital muscle disease. Molecular dynamics simulations suggested that the pathomechanisms by which the variants cause muscle disease include disruption of the binding sites for Ca2+ and for troponin I. In line with these findings, physiological studies in myofibers isolated from patients' biopsies revealed a markedly reduced force response of the sarcomeres to [Ca2+]. This pathomechanism was further confirmed in experiments in which contractile dysfunction was evoked by replacing TnC in myofibers from healthy control subjects with recombinant, mutant TnC. Conversely, the contractile dysfunction of myofibers from patients was repaired by replacing endogenous, mutant TnC with recombinant, wild-type TnC. Finally, we tested the therapeutic potential of the fast skeletal muscle troponin activator tirasemtiv in patients' myofibers and showed that the contractile dysfunction was repaired. Thus, our data reveal that pathogenic variants in TNNC2 cause congenital muscle disease, and they provide therapeutic angles to repair muscle contractility.

Long QT Syndrome Type 2: Emerging Strategies for Correcting Class 2 KCNH2 (hERG) Mutations and Identifying New Patients.

Significant advances in our understanding of the molecular mechanisms that cause congenital long QT syndrome (LQTS) have been made. A wide variety of experimental approaches, including heterologous expression of mutant ion channel proteins and the use of inducible pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) from LQTS patients offer insights into etiology and new therapeutic strategies. This review briefly discusses the major molecular mechanisms underlying LQTS type 2 (LQT2), which is caused by loss-of-function (LOF) mutations in the KCNH2 gene (also known as the human ether-à-go-go-related gene or hERG). Almost half of suspected LQT2-causing mutations are missense mutations, and functional studies suggest that about 90% of these mutations disrupt the intracellular transport, or trafficking, of the KCNH2-encoded Kv11.1 channel protein to the cell surface membrane. In this review, we discuss emerging strategies that improve the trafficking and functional expression of trafficking-deficient LQT2 Kv11.1 channel proteins to the cell surface membrane and how new insights into the structure of the Kv11.1 channel protein will lead to computational approaches that identify which KCNH2 missense variants confer a high-risk for LQT2.

The Role of a Crystallographically Unresolved Cytoplasmic Loop in Stabilizing the Bacterial Membrane Insertase YidC2.

YidC, a bacterial member of the YidC/Alb3/Oxa1 insertase family, mediates membrane protein assembly and insertion. Cytoplasmic loops are known to have functional significance in membrane proteins such as YidC. Employing microsecond-level molecular dynamics (MD) simulations, we show that the crystallographically unresolved C2 loop plays a crucial role in the structural dynamics of Bacillus halodurans YidC2. We have modeled the C2 loop and used all- atom MD simulations to investigate the structural dynamics of YidC2 in its apo form, both with and without the C2 loop. The C2 loop was found to stabilize the entire protein and particularly the C1 region. C2 was also found to stabilize the alpha-helical character of the C-terminal region. Interestingly, the highly polar or charged lipid head groups of the simulated membranes were found to interact with and stabilize the C2 loop. These findings demonstrate that the crystallographically unresolved loops of membrane proteins could be important for the stabilization of the protein despite the apparent lack of structure, which could be due to the absence of the relevant lipids to stabilize them in crystallographic conditions.

Lipid-Dependent Alternating Access Mechanism of a Bacterial Multidrug ABC Exporter.

By undergoing conformational changes, active membrane transporters alternate between an inward-facing (IF) and an outward-facing (OF) state to transport their substrates across cellular membrane. The conformational landscape of membrane transporters, however, could be influenced by their environment, and the dependence of the alternating access mechanism on the lipid composition has not been understood at the molecular level. We have performed an extensive set of microsecond-level all-atom molecular dynamics (MD) simulations on bacterial ATP binding cassette (ABC) exporter Sav1866 in six different phosphocholine (PC) and phosphoethanolamine (PE) lipid membrane environments. This study mainly focuses on the energetically downhill OF-to-IF conformational transition of Sav1866 upon the ATP hydrolysis. We observe that the transporter undergoes large-scale conformational changes in the PE environment, particularly in the POPE lipids, resulting in an IF-occluded conformation, a transition that does not occur when the transporter is embedded in any of the PC lipid bilayers. We propose that the PE lipids facilitate the closing of the protein on the periplasmic side due to their highly polar headgroups that mediate the interaction of the two transmembrane (TM) bundles by a network of lipid-lipid and lipid-protein hydrogen bonds. POPE lipids in particular facilitate the closure of periplasmic gate by promoting a hinge formation in TM helices and an interbundle salt bridge formation. This study explains how the alternating access mechanism and the flippase activity in ABC exporters could be lipid-dependent.

What Can and Cannot Be Learned from Molecular Dynamics Simulations of Bacterial Proton-Coupled Oligopeptide Transporter GkPOT?

We have performed an extensive set of all-atom molecular dynamics (MD) simulations of a bacterial proton-coupled oligopeptide transporter (POT) in an explicit membrane environment. We have characterized both the local and global conformational dynamics of the transporter upon the proton and/or substrate binding, within a statistical framework. Our results reveal a clearly distinct behavior for local conformational dynamics in the absence and presence of the proton at the putative proton binding residue E310. Particularly, we find that the substrate binding conformation is drastically different in the two conditions, where the substrate binds to the protein in a lateral/vertical manner, in the presence/absence of the proton. We do not observe any statistically significant distinctive behavior in terms of the global conformational changes in different simulation conditions, within the time scales of our simulations. Our extensive simulations and analyses call into question the implicit assumption of many MD studies that local conformational changes observed in short simulations could provide clues to the global conformational changes that occur on much longer time scales. The linear regression analysis of quantities associated with the global conformational fluctuations, however, provides an indication of a mechanism involving the concerted motion of the transmembrane helices, consistent with the rocker-switch mechanism.

A review of monoamine transporter-ligand interactions.

Transporters of the monoamines serotonin, dopamine, and norepinephrine are plasma membrane proteins belonging to the neurotransmitter sodium symporter family (NSS). Monoamine transporters (MATs) by facilitating reuptake of neurotransmitters from the synapse into the presynaptic nerve terminal, regulate neurotransmitter chemical signaling and maintain homeostasis. MATs are targets for several psychostimulants such as cocaine and amphetamine; and also for drugs treating several psychiatric disorders such as depression, attention deficit hyperactivity disorder, Parkinson's disease, and schizophrenia. Since, currently available treatment has several limitations and side effects, novel treatment is highly desired. Efforts to develop better treatment have been hampered by the lack of crystal structures for MATs. However, leucine transporter (LeuTAa), a bacterial protein from Aquifex aeolicus, belonging to the same NSS family as MATs has recently been crystallized. LeuTAa is used as a template to develop homology models of MATs, which facilitates understanding of the structure, function and pharmacology of MATs. Experimental methods for drug discovery demand a significant amount of time, effort and money. Efficient utilization of computational techniques hastens the process of drug discovery and also significantly reduces the cost. Assessing the binding affinity of drugs to the receptors is a key aspect of drug design. Free energy calculations compliment the experiment by quantitatively assessing the affinity of ligands to receptors. These methods are highly beneficial in the lead identification and optimization stages of rational drug design. We review the currently available free energy methods to treat protein-ligand interactions along with several free energy studies performed on MATs.

New design strategies for antidepressant drugs.

INTRODUCTION: In spite of research efforts spanning six decades, the most prominent antidepressant drugs to date still carry several adverse effects, often serious enough to warrant discontinuation of the drug. Molecular mechanisms of depression are now better understood such that some of the specific receptors responsible can be targeted for activation or inhibition. This advance, coupled with the recent availability of crystal structures of relevant drug targets or their homologs, has opened the door for new antidepressant therapeutic compounds. AREAS COVERED: The authors review the evolution of monoamine-based antidepressant drugs, up to the selective serotonin reuptake inhibitors (SSRIs). The authors discuss classic and contemporary antidepressant drug design strategies, with a focus on virtual screening and fragment-based drug design methods. Furthermore, they discuss the recent advancements in the understanding of the serotonin transporter (SERT) structure/function relationship in the context of recognition of SSRIs and outline a strategy for the use of computational approaches in producing new SSRI lead compounds. EXPERT OPINION: The authors suggest that given the long-awaited availability of credible three-dimensional structures for the SERT and related monoamine transporter proteins, cutting-edge computational methods should be the linchpin of future drug discovery efforts regarding monoamine-based antidepressant lead compounds. Because these transporter inhibitors cause a ubiquitous increase in extraneuronal neurotransmitter levels leading to side and adverse therapeutic effects, the drug discovery should extend to appropriate manipulation of the 'downstream' receptors affected by the neurotransmitter boost. Efficient use of new computational strategies will accelerate the drug discovery process and reduce its economic burden.